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We give a criterion to differentiate between dissipative and diffusive quantum operations. It is based on 
the classical idea that dissipative processes contract volumes in phase space. We define a quantity that can be 
regarded as "quantum phase space contraction rate" and which is related to a fundamental property of quantum 
channels: non-unitality. We relate it to other properties of the channel and also show a simple example of 
dissipative noise composed with a chaotic map. The emergence of attaractor-like structures is displayed. 
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I. INTRODUCTION 

Quantum dissipative processes in the context of quantum 
optics d, of superradiance Ql, of cooling mechanisms in 
ion traps |3] and open nanostructures |4] have been widely 
studied using a master equation approach. The corresponding 
master equation for each case is derived by modeling a micro- 
scopic interaction between the system (according to the case: 
an oscillator, a large spin, etc.) and a heat bath representing 
the environment. Under standard approximations (e.g. Born- 
Markov, random wave) (jl one ends with a Lindblad master 
equation, which automatically insures the preservation of pos- 
itivity |5J,|g]. For general dissipative processes the operators 
in the master equation are not normal. 

In this article dissipation is understood in analogy with 
classical dynamics, that is, we call dissipative any process in 
which phase space volume is not preserved. We have to keep 
in mind that this definition is different from the concept of 
'dissipativity' which is often found in the literature. For ex- 
ample, in 1 5] Lindblad calls dissipative maps the generators of 
completely positive semigroups, that is, dissipativity is linked 
to irreversibility. For Lidar et al.|7|] the condition of dissi- 
pativity is related to strictly purity decreasing, and shows its 
close relation to unitality. The concept of phase space con- 
traction is well defined in a classical context as the divergence 
of the velocity field. In quantum mechanics this divergence 
can be related to the sum of the commutators of the Lindblad 
generators in the classical limit. |8] 

The evolution of the density matrix of the system between 
two given times, that is, the transformation that takes p(f) to 
p(f + 8t) can be obtained by integration of the master equa- 
tion. This leads to a superoperator which can be written in an 
operator sum or Kraus form |9]. The detailed procedure de- 
pends on each particular process and can be mathematically 
complicated. In this work we follow an alternative approach 
to the description of dissipative quantum noise which consists 
in directly modeling the superoperator in its Kraus represen- 
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tation. The formalism of quantum operations to describe open 
systems, which is especially well adapted in the context of 
quantum information and quantum computin g llClJTlll . has 
also been used for open quantum maps | !^ [T3lfT4lll5LllrjlfT7ll . 
which are quantum maps in the usual sense with some non- 
unitary noise that can be considered the effect of an interaction 
with some environment. 

In 11511 the authors analyze several models for purely dif- 
fusive noises and show that generalization of simple chan- 
nels known in quantum information theory (depolarizing and 
phase damping channels) can be written as an incoherent sum 
of translations in phase space. Since the translation operators 
are unitary the preservation of the trace implies that the noise 
superoperator is unital (i.e. identity preserving). In this pa- 
per we focus on non-unital channels in a finite dimensional 
Hilbert space and show that their action can be interpreted as 
a contraction of the phase space volume leading to dissipative 
dynamics. 

This paper is organized as follows. In Sec.|II|we briefly go 
through the formalism of the Lindblad master equation and the 
superoperator approach. A parameter measuring phase space 
contraction in the quantum operations formalism is introduced 
in Sec. [II]] Some examples of quantum diffusive and dissipa- 
tive processes are given inSec HVl In Sec.|V]we give a simple 
dissipative channel as an example and show the effects when 
composed with a unitary map in Sec. IVII 



II. MASTER EQUATION AND QUANTUM DYNAMICAL 
SEMIGROUPS 

Open systems involve non-unitary evolution of states into 
probabilistic ensembles of states, or mixed states. The de- 
scription of mixed states is given by the density operator p . 
Usually the dynamics of open quantum systems is described 
by a Hamiltonian 

H = H s +H r + Hi (1) 

consisting of the dynamics of the system (s) and the reser- 
voir (r) and an interaction term. If the Markov approximation 
holds for the reservoir, then the problem can be reduced to a 
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master equation for the density matrix of the system, 



(2) 



which in general can be written (in the interaction picture) in 
GKS form |6] and in turn can be simplified into the so-called 
Lindblad |5] form 

| = ^I{^p4]+[ap4]} (3) 

where the Lindblad operators L, are system operators which 
are in principle arbitrary and which characterize the noise. 
The formal solution of the master equation is 



p(At) 



„Af_Sf 



P(0) 



(4) 



so Jz? is the generator of a one-parameter family of com- 
pletely positive (CP), trace preserving (TP) linear maps S 
called quantum dynamical semigroup\5, 6]. These operators 
acting on the space of density matrices (the subspace of pos- 
itive semidefinite operators include the density matrices) are 
given such names as superoperators, quantum operations or 
channels depending on the context in which they appear. We 
use any of those names indistinctively. Moreover in this work 
we consider only discrete time CP maps which can (but not 
necessarily do) arise from an integration of a Lindblad equa- 
tion through time Af . 

It is well known fiol Hill that a completely positive super- 
operator can be written in the operator sum or Kraus 1 9] form. 
That is, there exists a set of operators such that 



S(p) = £M M pM^ : 



(5) 



where p. < N 2 , with the dimension of the state space, the 
map S is TP if 



(6) 



In the analysis that follows pure states are vectors in finite 
dimensional Hilbert space ,J#n (of dimension N) and therefore 
density operators are represented as complex A x N matrices. 
A finite Hilbert space of dimension N is naturally associated 
to a classical phase space of finite volume (that we normal- 
ize to unity). The dimension is then given by the semiclas- 
sical rule N = l/h The volume occupied by a state p is well 
represented by the purity Tr[p 2 ]. The simplest examples are 
the sphere (for angular momentum applications) and the torus 
(for quantum maps). The association of quantum density ma- 
trices to distributions in phase space can be done in several 
well known ways via the Weyl-Wigner, Husimi or Kirkwood 
representations. 



III. QUANTUM PHASE SPACE CONTRACTION 

The standard definition of dissipation in classical mechan- 
ics is related to the contraction of phase space volumes given 



by the divergence of the drift vector entering the Fokker Plank 
equation | jjj. 

At the level of the quantum master equation it is well un- 
derstood how the properties of the Lindblad operators deter- 
mine whether dissipation will be present in the correspond- 
ing classical dynamics Hil Il9ll . In (2(1 Strunz and Percival 
derive a Fokker Planck equation by taking the semiclassical 
limit of a quantum master equation of the type and show 
how non-Hermitian Lindblad operators lead to classical dissi- 
pation (and quantum fluctuations), while Hermitian ones de- 
scribe diffusion (these could also have a contribution to dissi- 
pation, but in usual physical applications they do not). They 
obtain an expression to lowest order in h of the phase space 
divergence of the vector drift A in terms of the Poisson brack- 
ets of the Wigner-Moyal transforms Lj(q,p) of the Lindblad 
operators 



divA = -iY,{Li{q,p),Lj(q,p)}. 



(7) 



Given the correspondence between Poisson brackets and com- 
mutators it is clear from this equation that dissipative pro- 
cesses (divA 7^ 0) correspond to non normal Lindblad oper- 
ators (i.e. LiL] = l]L{) in the master equation 0, that is to 
non-unital processes which by definition satisfy «5f (/) ^ or, 
at the level of superoperators S(I) ^ I. 

Dissipation is then described by non-unital quantum opera- 
tions. On this basis we now introduce a dissipation parameter, 
which measures the non unitality of a given quantum channel, 



Tr 



(S( Pl )- P iY 



Tr[(p/ 



= NTr[F 2 }, 



where p/ = I/N is the maximally mixed state and 

r=S(p / )-p / =i£[M M ,Mt] 



(8) 



(9) 



is a traceless, Hermitian operator which gives a measure of the 
non-normality of the Kraus operators . 

Dissipation leads to a concentration of probability in some 
states (with a consequent reduction in others). In this sense, 
a non-vanishing 77 implies that a contraction of the available 
phase space volume with respect to the uniform distribution 
has been achieved. It is this contraction that is measured by 
the parameter 77 . 

This operator has significance in a phase space representa- 
tion. For example, its Husimi function 



(ziriz) 



(10) 



where \z) is the coherent state centered at z — (q,p) (and z* 
implies complex conjugate), gives a local description of the 
dissipation process. However, this is not the only possibility 
as other operator bases may be used to provide different pic- 
tures. For example, for quantum information applications the 
Pauli 1 10] basis plays a prominent role 1 30]. 
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A more precise connection between the parameter tj and a 
classical contraction rate can be made by writing Eq. (|8j in 
terms of the variation of the purity 



which leads to 



p„(p)=Tr[p 2 ]=Tr[(S"(p)) 2 



(11) 



at "time" n. Expanding the square in Eq. (|8j it follows that 
Tr[S( P/ ) 2 -( P/ ) 2 ] _pi( P i)-po(pi) 



V = 



Tr[(p/) 2 



Po(Pi) 



(12) 



This quantity can also be related to the Lindblad operators 
by taking At in Eq. (0} small enough so that a quadratic ex- 
pansion of the purity 



Pt(pi) = PoiPl) +po(pi)At + l -po{pi)At 2 



(13) 



is valid. Since the first derivative of the purity for p/ vanishes, 
we find (for Af = 1) 



Pi{Pi)-po(p, 



Po(Pi 



- = il Tr (Wj] [Lj,L]] 



In the classical limit, from Eqs. @ and (|8j we get: 

(div A) 2 dqdp 



n 



(14) 



(15) 



and thus tj can be understood as a global contraction rate. 

There is another way to interpret the parameter tj. The map 
S acts linearly on the space of operators. Therefore it has asso- 
ciated a matrix representation. If we consider an orthonormal 
basis of N 2 operators A,-, such that Ao = I/VN and the rest 
are traceless operators, the matrix representation of S takes 
the form 



" 1 


V2 








(16) 


Vl 


M 





where 



[% = Tr[AjS(A;)]. 



(17) 



The matrix M is (N 2 — 1) x (N 2 — 1) and Vi and V2 are col- 
umn and row vectors respectively (we use the square brackets 
to identify the representation of a superoperator or an opera- 
tor in this basis). This is the so-called affine representation. 
The trace preserving condition is met if V2 = while v i = 
implies that the map is unital. So if V] ^ in Eq. (II 61 . the 
contraction parameter is 



n 



[p/] f [p/] 



(18) 



Now in this basis 



[Pl} = 



fl/VN\ 




V o J 



(19) 



77 = vjvi 



(20) 



and the contraction information is contained in the first col- 
umn of [S] . In a way tj measures the non-symmetric form the 
matrix takes in the affine representation. 

It is worth remarking that while the left invariant eigenstate 
of [S] is [p/] the right invariant eigenstate is something differ- 
ent which in the case of the noise composed with a unitary, 
possibly chaotic, map can take a very involved form resem- 
bling the fractal structure of a strange attractor. An example 
of this is given in Sec. lVIl 



IV. EXAMPLES 

In this section we analyze some examples of purely dif- 
fusive and dissipative noise models found in the literature to 
illustrate the classification introduced above and compute the 
corresponding contraction parameter tj. 

Our definition does not make any explicit reference to 
the classical phase space, it only refers to finite dimensional 
Hilbert space. In what follows we choose to use the torus 
which can be represented as a square of unit side, with peri- 
odic boundary conditions. Therefore the superoperator S acts 
on the space of complex N x N matrices. 



A. Random unitary processes 

Trace preserving, unital processes are also called bistochas- 
tic. We identify bistochastic quantum operations, with non 
trivial Kraus representation, with pure diffusion (possibly with 
drift) while we associate dissipation with non-unital maps, 
characterized by a contraction parameter tj 7^ 0. 

A typical example of bistochastic map is a random sum of 
unitary operations 



s(p) = Y,c^ P uj l , 

with Ufi unitary and the trace preserving condition 

Y< C H = 1' C M ^ °> 



(21) 



(22) 



also known as random unitary proces s (RUP) . Diffusive noise 
in the form of a RUP was studied in fl3lfl4lfl7il . as a Gaus- 
sian sum of (normalized) translations in phase space. When 
composed with a unitary map the whole noisy map can be in- 
terpreted as a coarse graining ll4ll7ll of the original map. 

Of course there can be other examples of bistochastic (or 
purely diffusive) channels which need not be RUP's. The 
only necessary condition (which follows from Eq. (|5Jl) is that 
the Kraus operators be normal. In the context of quantum in- 
formation noise is characterized by operations that can either 
take place on single or multiple qubits, and consist in com- 
binations of (tensor products of) Pauli matrices. While some 
well known noises can be re-expressed as RUP's 1 1511. it is not 
the most general situation 1I2II1 . 
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B. Sloppy Bakers 

Recently a model of an irreversible quantum baker map was 
presented by Lozinski et al. fifill . It consists of two steps: the 
usual quantization !^ l23ll of the baker map which is unitary, 
and a projective measurement with a controlled translation in 
momentum. That is the bottom half of the torus is left un- 
touched while the top half is translated an amount A/2 The su- 
peroperator is non-unital except for the case where A = 0. The 
invariant state is a uniform distribution on a reduced Hilbert 
space of dim = N(l — A). Thus the original invariant state of 
the unitary baker's map is reduced by a strip of area A (see 
Fig. 2 in CI). 

We do not reproduce the exact expression for the superop- 
erator but the parameter 77 for this map can be shown to be 

77 = A. (23) 

As expected 77 depends directly on the model's contraction 
parameter A in a simple way. 
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FIG. 1 : Sechematic representation of the Husimi function T(z, z* ) for 
the noise used in 1 16]. Dissipation takes place where T(z,z*) 7^ 0. 

The operator A defined in Eq. (|9} for this model is 

r= l(y-NA/2 Dt yNA/2_ D ^ (24) 

where V NA / 2 is a translation size A/2 in momentum and D t 
is the projection onto the upper half of the torus. For clar- 
ity in FIG. [D instead of the Weyl representation of F, the 
corresponding Husimi T(z,z*) is shown. The regions where 
r(z) 7^ are the regions where dissipation takes place. 



C. Generalized amplitude damping noises 

A large variety of dissipative noise models which are de- 
rived from a microscopic Hamiltonian for a system (an oscil- 
lator 1 24], a rotor |25], a large spin |2, 26], etc) in contact 



with a bath in thermal equilibrium can be interpreted as gen- 
eralized amplitude damping models and modeled by the fol- 
lowing Kraus superoperator: 

S(P)=I>P A ! (25) 

k 

where 

A k = Y,c k n \n)(n + k\ (26) 
is a combination of the transition operators 

Pij = m\- (27) 

Although the basis states |/) and the form of the coefficients 
depend on each particular problem, a common feature to 
all these models is that when S acts on the basis of the Py the 
skewness a = (i — j) is conserved. With the notation Pij = Pf 
the evolution of the P, a in each subspace labeled by a is given 
by 

S/f=L</^ (28) 

/ 

withm" =c i r l c i l Z l *. 

The preservation of the trace implies that Y^=o m ?/ = 
1 while the unitality condition corresponds in addition to 
Y^Jq 1 wtfi — 1 ■ It is then clear that in order to analyze the 
unitality of a noise of this type it is sufficient to study the 
properties of the (N x A^)- matrix m%, unitality corresponding 
to the bistochasticity of this matrix. 

It is immediate to see that if the matrix m?, is either sym- 
metric, or symmetric with respect to the diagonal i + j = N— 1 
stochasticity implies bistochasticity. 

A symmetric matrix corresponds to having a purely diffu- 
sive reservoir, i.e., a thermal bath at infinite temperature, ex- 
changing quanta with the system in both directions at equal 
rate. 

In the second case the matrix elements only depend on 

(i — j), implying that in Eq. (I26> the coefficients = c(k). 
Expanding the transition operators in the basis of the transla- 
tions on the torus Tj^.p) (as defined by Schwinger|27]), it is 
easy to see that the superoperator S is diagonal in this basis 

S P = pL c (-*)W7jU) (29) 

that is, it can be written as an incoherent sum of translations 
in phase space. 

If we compute 77 for this type of models we get: 

T7=L(I>,/) 2 -l- (30) 

/ i 

For the dissipative map studied by Dittrich and Graham in 
ll25ll (adapted to a finite dimensional Hilbert space), this gives 
in the limit of zero temperature (to second order in 7) 

Tj=y 2 (iV-l) ) (31) 

where 7 is the friction parameter (differential rate of loss of 
action) of the corresponding classical map (see Eq. (2.4) in 
10). 
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V. SIMPLE DISSIPATION CHANNEL 

We propose a family of non-unital noise channels to serve 
as simplified models of dissipative processes like the one de- 
scribed in Sec. lIV CI Let 



(32) 



D e (p) = (l-e)p + e£p iJ PijpPf j , 



with Pij as in Eq. J27t and py real and positive. Equation (|3^ 
is in Kraus form and thus D £ is CP. The TP condition is 



(33) 



where T means 'transpose'. Therefore the matrix of coeffi- 



\ 




FIG. 2: Graphical representation of the structure of the matrix [py] 
with N = 32 and a = 0.5. The color scheme is as follows: white 
means null value and a black square equals one element of value 
equal to 1. 
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FIG. 3: (Color online) The contracting parameter T) as a function of 
a for the dissipation model of Eq. 1341 for two different values of e: 
(*) e = 0.5, (A) e = 0.75. 



cients pij should be stochastic. If it is doubly stochastic then 
the noise is unital. This is a model which is diagonal in the pj 
representation and decoherence appears as a reduction by a 
factor (1 — e) of the non-diagonal terms of the density matrix. 



We can create dissipative models just by using non- 
symmetric stochastic matrix of coefficients pij. Consider the 
family 



D £a (p) = (l-e)p + e 



N/2 
i=-N/2 



iP P \ai], 



(34) 



with a E [0, 1) and the index [ai] is the integer part of ai, as a 
function of the parameter a and the negative values are taken 
mod(Af). In FIG.|2]there is an illustration of the structure of 
the matrix [py] of Eq. i32\ for the case a = 0.5, N = 32. 

For all values of a, the first term attenuates uniformly the 
off-diagonal elements of the density matrix, introducing de- 
coherence. The second term which acts on the probabilities 
(diagonal elements) takes the system to lower states. The per- 
mitted transitions of the diagonal elements are determined by 
the parameter a but its dependence on the parameter e, which 
accounts for the coupling strength of the system to the envi- 
ronment, becomes important when the contraction is consid- 
ered. 

If Pij are transitions in the computational basis of a state 
of say k = \og 2 N qubits, then the noise can be interpreted 
as follows. With probability (1 — e) it leaves the initial state 
untouched while with probability e it induces errors in the 
form of transitions which depend on the parameter a. For 
a = 1 the superoperator is unital and corresponds to a phase 
damping channel for k = \og 2 N qubits . If a < 1 this very 
simple model captures some important features of the am- 
plitude damping noises (an example of which is described in 
Sec.lE - 



On the other hand, if we take Py to be transitions 
in momentum and since the invariant state (for all a ^ 0) is 
p* = Poo = \p = 0){p = 0| then the noise has the effect of a 
friction. As is shown later, in phase space representation the 
parameter a is related to the region over which dissipation 
acts. 

The simple form of Eq. i34l allows to compute the contrac- 
tion parameter easily as 

77 = NTv[{D £a (I/N)-I/N) 2 ] 

-2 

-Tr /-2£^ cei])W +£/ , r, 

' ij 



N 



N 



-Tr 



El 



[aj\ai)\ z -N 



(35) 



where Y 1 iP[ai].[ai] = N an d f° r simplicity we drop the [ ] sym- 
bol inside the kets and bras. Now (aj\ai) = 0, 1 so the square 
can be dropped and in order to calculate the sum we approxi- 
mate by continuous variables as 

jj 8{a{xi-Xj))~- \l 8{xi- X j) 



an 



a 



a ij a 



Therefore we get 



T] = e 



1-a 



a 



(36) 



(37) 
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FIG. 4: Plot of the Husimi function F(z,Z*) for N = 64 (in gray scale, black is maximum and white minimum) and different values of a.Lighter 
regions indicate local contractive regions. 



In figure FIG.|3lve plot rj(a,e) for two different values of 
e. We see that the computed points fit exactly the analytic 
expression. Moreover the saturation value, 



77(a = l/N,e) = e 2 {N-l). 



(38) 



for small a can be understood as follows. Since the dimension 
of Hilbert space is N then the integer part [a i] (i = 0,...N—l) 
for any value of a < 1 /N is zero. We notice that the depen- 
dence on the coupling parameter e is the same as the one given 
in Eq. Ql> for the dissipative map of l25ll . if one identifies it 
with the friction parameter y. In addition, it can be easily seen 
by taking the mean values of p and q that this operation im- 
plies the following dissipative map in the classical limit 



(l-e)p 



(39) 



where the bar indicates mean values. 

For this model the operator F defined in Eq. l|9} is 



N/2 

' E I 

i=-N/2 



Pii) 



(40) 



So the region of dissipation (determined by a negative value 
of A) has area equal to 1 — a, and occupies the central region 
of the unit square. In figure FIG. the Husimi function of 
the operator F is represented taking different values of a (and 
£ = 0.5). The light region represents the area where contrac- 
tion takes place and corresponds to a negative value of F. For 
a ~ the contraction is uniform over (almost) all phase space, 
except over the state |0) of momentum. 



VI. COMPOSITION WITH A UNITARY PROCESS 

Following recent works 0H[HE[l7l El IU we study 
the effect of the dissipative noise channel described in Sec. |V] 
when composed with a unitary map. As an example we take 
the quantum version of the standard map on the torus. We 
suppose that to a good approximation the whole noisy propa- 
gation takes place in two steps 



where U(p) = UpU^ is the unitary step. The two-step scheme 
can can be used, in the master equation if the Hamiltonian 
part commutes, to some desired order in h with the non- 
Hamiltonian part. In the quantum operation formalism these 
scheme is suitable in the case for example where unitary evo- 
lution takes place in so short times that the noise is negligible 
(e.g. the micro-maser, a billiard where the interaction with 
the walls is very short and the evolution inside is dissipative). 
The unitary map chosen is the quantum version of the standard 




FIG. 5: (Color online) Left: histogram representation of the attractor 
for the classical standard map (Eq. J42j with k = 0.065, 5 = 0.6) with 
dissipation. Right: Husimi representation for invariant state for the 
quantum standard map with the dissipative noise described in Secfvl 
(a~0,£ = 0.4 



map on the torus 



q + p 

8p — 2nksm(2%q) 



(42) 
(43) 



S(p) = D ea (U(p)) 



(41) 



where the factor 8 < 1 on the momentum acts as friction. 

In FIG.|5]the classical and quantum invariant states are plot- 
ted. The classical attractor on the left is represented as density 
in phase space. On the right is the Husimi distribution of the 
invariant state for the quantum version of the map followed by 
the noise of Eq. f34Tl for a ~ 1 /N (N = 64). Since the dissipa- 
tion terms in Eqs. J39i and (I42t are the same when 8 = (1 — e), 
as expected the quantum invariant state exhibits the structure 
of the classical attractor. 
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FIG. 6: Husimi representation (same color scheme as FIG.|4j of the invariant state for the standard map with dissipation for different values 
of a and £ and N = 64. The difference between a ~ (dissipation) and a ~ 1 (diffusion) can be observed. (Axes and units are the same as in 
FIGS.0and FIG.0) 



The shadow on the lower and upper edge of the torus can 
be explained from the definition of F, Eq. ( I40t . for this noise 
(corresponding to the rightmost image in FIG.|4). The same 
argument can be used to explain FIG.|6] As a grows the re- 
gion of dissipation becomes smaller as well as the contract- 
ing parameter tj . The noise becomes increasingly similar to 
a unital operation, at least on the upper and lower bands of 
width a/2. So over the black shaded regions of FIG. |6] the 
noise acts like a generalized phase damping channel fl5ll . 
where the preferred basis are the momentum projectors Pa 
with i = {0, . . . ,aN/2}(J {(1 -a/2)N- ,N—l}, while 
on the lighter region, a dissipation of the type of Eq. d39l acts 
and the attractor is uncovered. 



blad operators to the vector drift in its Fokker-Planck limit, 
we defined a parameter that measures non-unitality and char- 
acterizes dissipative quantum operations. As an example we 
proposed a noise channel that displays in a simple way the 
essential features of decoherent and dissipative processes. 

The non-unitality of the superoperator is related to a trace- 
less Hermitian operator whose phase space distribution gives 
a local image of the dissipation process. This operator is inde- 
pendent on the Kraus representation and can give a useful in- 
sight into the dissipative properties of the superoperator when 
a phase space description or a semiclassical limit is not avail- 
able, like for general noise channels in quantum information. 



VII. CONCLUSIONS 

The fundamental difference between unital and non-unital 
processes was explored. In analogy to the classical limit of the 
master equation, which relates the commutator of the Lind- 
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